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Abstract. In addition to the usual linear gravitational waves in transverse-traceless coordinates, 
higher-order gravitational field equations, as well as their corresponding solutions, are explicitly 
obtained. It is found that higher-order waves do not represent corrections to the hrst-order wave. 
In contrast, all higher than second-order solutions do represent corrections to the second-order 
wave, a property that makes the hrst-order gravitational wave to stand apart from higher-order 
waves. Furthermore, although the hrst-order solution is transverse and traceless, all higher-order 
solutions are not. As a consequence, the whole solution is neither transverse nor traceless, a 
result that could eventually have important consequences for quantum gravity, and in particular 
for the dehnition of graviton itself. Some additional properties and features of these higher-order 
gravitational waves are explored and discussed. 


1 Introduction 

Although there seems to be a general agreement that the transport of energy by gravitational 
waves is a non-linear phenomenon [l|, due to the tiny amount of energy transported, it is 
usual to assume that the linear approximation is able to accurately describe the propagation 
of gravitational waves j^. For this reason, when studying gravitational waves far away from 
the sources—that is, in the so-called wave zone—only the linear approximation is considered. 
Denoting by the first-order term of the post-Minkowskian expansion, and choosing the 
class of harmonic coordinates, in which satisfies the condition 

duK^)u - \duh-L) = 0 , ( 1 ) 

with /i(i) = the first-order vacuum field equation assumes the form 

= 0 . ( 2 ) 

Gravitational waves are usually interpreted as a solution to this set of equations. 

On the other hand, higher-order gravitational waves are in general studied through the use 
of an induction process, which amounts to consider non-linear iteration of the field equations. 
Such approach yields, for any n > 2, the set of higher-order equations 

^U^{n)u ~ = 0 (3) 

■ ■ ■ ^hn-1)]- W 

Even though the reneral mathematical properties of such system have been extensively explored 
in the literature [^, a thorough study of the explicit form of the field equations for n > 2, as 
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well as their corresponding solutions, is still lacking. With the purpose of getting a clearer and 
more precise view of higher-order gravitational waves, we are going in this paper to explicitly 
obtain a few higher-order held equations along with the corresponding coordinate conditions. 
Then, by solving these equations, some features and properties of higher-order gravitational 
waves will be studied and discussed. 


2 Gravitational waves 


2.1 Perturbative expansion 

Considering that the sources of gravitational waves are at enormous distances from Earth, 
it is sensible to assume that the amplitude of a gravitational wave when reaching a detector 
on Earth will be very small. This allows the use of a perturbative analysis in which the 
gravitational held variable—that is, the metric tensor—is expanded in the form 




(5) 


where e is a small parameter introduced to label the successive orders, and is the metric of 
the background Minkowski spacetime, where the gravitational wave is supposed to propagate. 
The expansion of the inverse metric is 

- e hH’; + + • • • . (6) 

At the hrst order in the metric perturbation, the Christoffel connection assumes the form 


^( i )/37 ~ 2^ ^ (<9/3^( i ) p 7 + d-yh(^i)pp dphi^i)jj^) 
whereas the hrst-order Riemann tensor is 

pa _ Q pa _ Q pa 

At the second order, the Christoffel connection is 

„ — In 
( 2)/37 2 7 

with the second-order Riemann tensor given by 


— n 'll ^ i dfiht2)p'y T d^h(2)pi3 dph(2)i3^ 2 ^(i) 


pa _ o pa _ 

^(2)pLVl3 {'2)p,0 ^p-*- (2)pi/ 

At the third order, the Christoffel connection is 


pa p7 _ pa p7 
(1)71^ (i)p.0 (1)7/3 WfJ-i' 


(7) 

( 8 ) 

(9) 

( 10 ) 


pa _ 1 ap 

^ (3)/37 — 2 'I 


3)P7 


T d'^^{3)p0 dph(3)fi'^ 2 ^( 1 ) ^^/3^(2)p7 T d'yh(2)p0 dp^(2)0^ 

T 2 (^( 2 ^ T A d^h(^i'jpis dph(^i'jjj^ 


and the third-order Riemann tensor given by 


pa _ Q pa _ 

And so on, to any higher order. 


1 pa p7 _ pa p7 
“T ^ (1)71/-^ {2)P0 ^ (1)7/3 (2)P!^ 


I pa p7 _ pa p7 

“T (2)71/-*- (I)p0 ^ (2)7/3^ (l)pi^' 


( 11 ) 

( 12 ) 
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2.2 First-order gravitational waves 

In the linear, or first-order approximation, the sourceless gravitational field equation is 

= 0 (13) 

with -R(i)^i/ the first-order Ricci tensor. In terms of the metric perturbation, it assumes the 
form 

^ = 0 (14) 

where the notation /i(i) = has been used. The class of harmonic coordinates is defined by 
the condition 

(7^"rv = o. (15) 

At the first order it becomes 

= 0- (16) 

In these coordinates, the field equation (11411 assumes the form 

□ h(V = 0 (17) 

with □ the flat-spacetime d’Alambertian operator. 

A monochromatic plane-wave solution to this equation is 

K^)v = (18) 

where is the symmetric polarization tensor. In this case, the wave vector is found to 
satisfy the dispersion relation 

kpkP = 0 (19) 

and the harmonic coordinate condition ()16p becomes 

= \kyh(^y ( 20 ) 

Within the class of harmonic coordinates, it is possible to choose a specific coordinate system. 

Once this is done, the coordinates become completely specified, and the remaining components 
of turn out to represent the physical degrees of freedom. A quite convenient choice is the 
so-called transverse-traceless coordinates (or gauge, in analogy to electromagnetism), in which 

/i(i) = 0 and hl^i)yU^o) = 0 (21) 

with an arbitrary, constant four-velocity. In this coordinate system, condition (I20p assumes 
the form 

kphl[^, = 0 . ( 22 ) 

Now, although the coordinate system has already been completely specified, we still have the 
freedom to choose different local frames e“ = e°‘pdxP. In particular, it is always possible to 
choose a specihc frame in which = ^o- frame, as can be seen from the second of 

the Eqs. (l2Tp . 

hP = 0 (23) 
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for all /i. Orienting the frame in such a way that the wave travels in the z direction, the wave 
vector assumes the form 

= (w/c, 0,0,w/c) , (24) 

and the physical components of the wave are found to be 

Km = -K)y and Km = hm""- (25) 

Linear waves satisfying these conditions are said to represent a plane gravitational wave in 
transverse-traceless coordinates. 

2.3 Second-order gravitational waves 

At the second order, the sourceless gravitational field equation acquires the form 

R{‘2)yv ~ 0 (26) 

with R(^ 2 )iiv the second-order Ricci tensor. In terms of the components of the metric perturba¬ 
tion, and using for K^^y the plane wave solution (fTSl) along with the constraints (fTOl) and (l20]l . 
the wave equation (j26p assumes the form 

- d^,{daK2)v - \dyh(2)) - dy{dah^2)y - kKK)) = “I exp[i2/cpX^] (27) 

where h( 2 ) = and 

4>(2) = (28) 

On the other hand, at this order the harmonic coordinate condition (|15p is 

= (29) 

Substituting Eqs. ([7]) and (l9|), it becomes 

dK(2)\ - \dxK) = -§$(2)fcAexp[i2A:px'']. (30) 

In these coordinates, and using for the plane wave solution (fT8]l . as well as the constraints 
(USD and (I20D, the second-order wave equation assumes the form 

n K)'^ ~ k^ky exp[i2kpX^]. (31) 

A monochromatic wave solution to this equation is 

= {Al'^^y + iBl^-^y) exp[i2A:pX^] (32) 

where 

= and B,%y = -^^k^^ky (33) 

are the second-order amplitudes, with Kg an arbitrary wave number four-vector. The ampli¬ 
tude of the hrst part of the solution satisfies 

^( 2 )g, — A( 2 ) = 4 ^( 2 ) and kpAj^^y = -g kyAg 2 ) (^4) 
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from where we see that, at this order, the solution is no longer traceless. On the other hand, 
the amplitude of the second part of the solution satisfies 


= 3 (^ 2 ) = 0 and = 0. 


(35) 


This part of the solution is thus traceless and transverse. 

We consider now a laboratory proper frame — endowed with a Cartesian coordinate sys¬ 
tem — from which the wave will be observed. If the wave is traveling in the z direction of the 
Cartesian system, the wave vector is that given by Eq. (IMl) . In this case, the amplitude 
turns out to be 




/I 0 0 0\ 

3 $( 2 ) 0 10 0 

16 0 0 1 0 

\0 0 0 1 / 


(36) 


with each component given by 


At — AX — A y — A ^ = A<i), , 

^(2)i ^(2)X -^(2)y ^(2)2 — 16 ^(2)- 


(37) 


On the other hand, choosing the arbitrary wave vector Kp in such a way that Kq = Ki = 
K 2 = 00 the amplitude turns out to be 


(-S{2)v) — 


/I 0 0 -1\ 
0 0 0 0 
8c 0 0 0 0 

\1 0 0 - 1 / 


(38) 


This part of the solution has two physical components, which are given by 

-^( 2)0 ~ ~k^(2)Z -^( 2)2 ~ ~B(2)Z^ ■ 


(39) 


2.4 Third-order gravitational waves 

Following the same steps of the previous section, we consider the third-order sourceless gravi¬ 
tational field equation 

R{ 3 )y.v = 0 (40) 

with R( 3 )py the third-order Ricci tensor. In terms of the components of the metric perturbation, 
and using the solutions already obtained for h^^^y and h^^^y, the field equation (|40]) assumes 
the form 

Oh^ 3 ^py - dp{dah^^y - \dy\ 3 )) - dy{dah“^^p - = - ^( 3 )kpky exp[i3kpxP] (41) 

where /i( 3 ) = hll^a and 

4*(3) = ^(i)o-^(i)A^(i)p- (42) 

‘This choice corresponds to assuming that the wave amplitude depends on the distance 2 from the 

source, but not on the time. 
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( 43 ) 


On the other hand, at this order the harmonic coordinate condition (|15jl is 

< h — n 

'/ (3)^1^ "'(1) ‘ (1)/^^^ ~ 

Substituting Eqs. (ED, m and EB, it becomes 

~ 5'9a^(2) = §^(3)A:a exp[i3/i;px^]. (44) 

In these coordinates, the third-order wave equation is found to be 

□ /i|^3)i/ = — 4'(3) k^ky exp[i3A:pX^]. (45) 

A monochromatic wave solution to this equation is 

exp[i3A:pX^] (46) 

where 

A 3 )v = ^^( 3 ) 5 (; and = (47) 

are the third-order amplitudes, with Kq an arbitrary wave number four-vector. There is a 
problem, though: it so happens that, as one can easily verify, $( 3 ) = 0 , and consequently 

= 0. (48) 

This means that the third-order solution does not contribute to the metric perturbation. As a 
matter of fact one can see that all odd-order terms in the metric perturbation vanish identically: 

^(5)!2 = ^(7)1^ = ^(9)1^ = • • • = 0. (49) 


2.5 Fourth-order gravitational waves 

At the fourth order, the sourceless gravitational field equation is 

Rwpv = 0 (50) 

with i?( 4 )pi/ the third-order Ricci tensor. In terms of the components of the metric perturbation, 
and using for and the solutions already obtained, the field equation ([50]) 

assumes the form 

- dp(da/i(" ,, - \dyh(^)) - di.(d„h(" p - = \ ^l^^k^ky expfid/cpx'’] (51) 

where /i( 4 ) = and $( 2 ) is given by Eq. (1281) . On the other hand, at this order the harmonic 
coordinate condition p5p is 

+ (< 2 ) + = 0 . (52) 

Substituting the Christoffel connections, it becomes 

dph{l^x - ^dxh(^4) = t^(2)^A exp[i4A:pa:^]. (53) 
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In these coordinates, the fourth-order wave equation assumes the form 


□ h'^i)u = - ^ ^( 2 ) exp[iAkpxP]. (54) 

A monochromatic wave solution to this equation is 

/i( 4 )i/ = exp[i4:kpxP] (55) 

where ^ 

= and = (56) 

are the fourth-order amplitudes, with Kq an arbitrary wave number four-vector. The amplitude 
of the first part of the solution satisfies 

= A( 4 ) = ^ and kpA^l-^i, = \ k^A^^^ (57) 

whilst the amplitude of the second part of the solution satisfies 

= . 6 ( 4 ) = 0 and kpB^^y = 0. (58) 


We consider now a laboratory proper frame — endowed with a Cartesian coordinate sys¬ 
tem — from which the wave will be observed. If the wave is traveling in the z direction of the 
Cartesian system, the wave vector is that given by Eq. (1241) . In this case, the amplitude A^^-^y 
turns out to be 


(-4(4)1/) - 4* ( 2 ) 


/I 0 0 0\ 
0 10 0 
0 0 10 
Vo 0 0 1/ 


(59) 


from where we see that this part of the solution has just one physical component: 

^{i)t — ^{i)x — ^(4)y — ^(4)Z — M *^( 2 )' 

On the other hand, choosing the arbitrary wave vector Kp in such a way that Kq = Ki = 
K 2 = 0 , the amplitude B^^^y turns out to be 


(-S(4)v) — 


/I 0 0 -1\ 

1 4>(22)2:a; 0 0 0 0 

32 c 0 0 0 0 

Vl 0 0 -1/ 


(61) 


From this expression we see that this part of the solution has two physical components: 


pO _ tdz 

■^( 4)0 ~ '^( 4)2 


and .B(° 4 )^ = -B( 4 )^°. 


(62) 


7 








2.6 Higher-order gravitational waves 


Based on the solutions found up to here, one can infer that the generic form of a higher-order 
solution is 

^Mi‘^nkpxp] (es) 

where 


and 


,^(2n)^) 


02n^(2) 


/I 0 0 
0 10 
0 0 1 
\0 0 0 


0\ 

0 

0 

1 / 



b2n^^ 


/I 0 0 
0 0 0 
0 0 0 
\1 0 0 



- 1 / 


(64) 


(65) 


with a 2 n and b 2 n numerical coefficients. Using Einstein’s equations and the harmonic coordi¬ 
nate conditions, given respectively by 


^{ 2 n)lJ.i^ — 0 and r^ 2 n)i^ — Oj 
these coefficients can be easily computed. The first few of them are 


06 = 
08 = - 
“10 = 


13 

4608 


be = - 


1 


83 


131072 

5023 

31457280 


^10 = 


1536 

65 

”147456 
163 

3145728 


( 66 ) 


(67) 


And so on to any higher order. Like the second and fourth-order waves, all higher-order 
solutions have three independent physical components. 


3 Higher-order curvature invariants 

In general relativity, it is possible to construct two invariants in terms of the components 
of the Riemann curvature tensorjll They are the Kretschmann (K) and the Pontryagin (P) 
invariants, defined respectively by 

K = and P = R^pp.^R^^P'^ 

Tn the general case there are actually three invariants, which degenerate in two for sourceless solutions to 
Einstein equation, which is our interest here. 














with pa the curvature Hodge dual. The first is a scalar whereas the 

second is a pseudo-scalar quantity. Now, as is well-known, for linear gravitational waves, both 
Kretschmann and Pontryagin invariants are found to vanish identically [^-@]: 

i^( 2 ) = = 0 and P( 2 ) = = 0. 

The question then arises: is this a specific property of linear gravitational waves, or it is 
a general property of waves of any order? Using the higher-order solutions obtained in the 
previous section, we are able to give an answer to this question. 

To begin with, using the first and second-order solutions, the third-order Kretschmann and 
Pontryaguin invariants are found to vanish identically, 

K,,) ^ = 0 ( 68 ) 

and 

^(3) = = 0. (69) 

Before considering the fourth-order invariant, it is important to remark that, although the 
third-order solution /i( 3 )^,y vanishes identically, the third-order Riemann tensor does not vanish 
because there are contributions to it coming from and h( 2 )py With this proviso, the 

fourth-order curvature invariants are then written as 

^(4) ^ + RwfSf.uRoJ^'^ (70) 

and 

p _ * po p fifiu . * po p fifiu . r>a * p /yi \ 

4 ^(4)— ^( 2 )fi^v^( 2 )oi ~r ^(3)a ■ 

Performing these computations, both fourth-order curvature invariants are also found to vanish 

K( 4 ) = 0 and P( 4 ) = 0. (72) 

The above pattern repeats for any higher order. For example, as already discussed, the fifth- 

order metric perturbation is found to vanish: /i( 5 )^j^ = 0. Using this result it is then possible 
to compute the fifth-order curvature invariants, which are found to vanish, 

K( 5 ) = 0 and P( 5 ) = 0. (73) 

Analogously, the sixth-order invariants can also be shown to vanish, 

K(e) = 0 and P(e) = 0. (74) 

We can then conclude that both Kretschmann and Pontryagin invariants vanish for gravita¬ 
tional waves of any order. 
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4 Discussion 


We have developed a systematic study of higher-order gravitational waves with the purpose 
of exploring their properties and the structure of the perturbation scheme. Considering that 
the first-order solution is computed in harmonic coordinates, and taking into account that 
such solution appears in all higher-order equations, the same class of coordinates must be used 
for dealing with higher-order gravitational waves. Furthermore, within this class, a specific 
coordinate system is then chosen, in terms of which the first-order perturbation becomes 
transverse and traceless. In this case, as is well-known, the harmonic coordinate condition is 
equivalent to the divergence-free condition ([T]). However, as we have explicitly verified, at orders 
higher than one the harmonic coordinate condition does not generalize to Q , but to expressions 
with non-vanishing right-hand side. For example, at the second-order the harmonic coordinate 
condition assumes the form (130 jl . The use of the correct harmonic coordinate condition is 
crucial for the mathematical consistency of the higher-order solutions. 

An intriguing result of these higher-order solutions is that the second-order solution does 
not represent corrections to the first-order wave, but rather shows up as a completely different 
solution. In contrast, all higher than second-order solutions do represent corrections to the 
second-order wave. The first-order gravitational wave, therefore, stands apart from higher- 
order waves. This picture bears a strong resemblance to the case of nonlinear surface waves in 
shallow water, where solitary waves can be obtained from the Navier-Stokes equation (for an 


inviscid fluid) through a perturbation scheme 10|]. At first order one obtains a linear wave- 


equation whose solution does not represent a physical wave; its unique role is to determine 
the dispersion relation of the system. At the second order, one obtains a nonlinear evolution 
equation—the so-called Korteweg-de Vries equation—whose solution represents the solitary 
wave observed in nature. Like in the gravitational case, the second-order solution does not 
represent a correction to the first-order solution, but shows up as a completely different wave. 
Higher-orders terms of the perturbative scheme, however, do yield corrections to the second- 
order wave. 

Another interesting feature of the higher-order perturbation analysis is that, except for the 
first order, all odd-order terms of the metric perturbation vanish identically. This feature may 
be related to the physical properties of the gravitational interaction. In fact, it is well-known 
that, for a spin-2 gau ge b oson, the particle-particle static interaction is always attractive in the 
linear approximation 1]|. On the other hand, if presented in the perturbative series, odd-order 
terms would contribute with a repulsive static interaction between particle-antiparticle. That 
the perturbative series of Einstein equation naturally exclude all odd-order terms of the expan¬ 
sion can then be considered an evidence that the static interaction between matter-antimatter 
is always attractive at any order. It should be remarked that, in spite of compelling arguments 
favouring the idea that the static particle-antiparticle gravitational interaction should also be 
attractive, this is as yet an open question in the sense that it has never been verified experi¬ 
mentally. We mention in passing that there are currently some experiments aiming exactly at 
probing the gravitation behaviour of antimatter. One example is the ALPHA experiment at 
CERN 


1 ^. 


A crucial output of the perturbation scheme is that, although the first-order solution is 
transverse and traceless, all higher-order solutions are not. In fact, the generic form of the 
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higher-order metric perturbation is 


^M^'^nkpxp] ( 75 ) 

with 

KqT^ 

^{2nY ~ ®2n‘^’(2) K ^(2n)^ ~ ^2n ^(2) 

We can thus clearly see that, while the i?-part is transverse and traceless, the A-part is neither 
transverse nor traceless: 


k„ = 02^^^(2) ki, 


^{2n)^ 


and 


A 


(2n)M 


= 4a2n^ 


( 2 )- 


(77) 


As a consequence, the whole solution, that is, the solution that includes higher-order con¬ 
tributions, is neither transverse nor traceless. This result may eventually have important 
consequences for quantum gravity, and in particular for the definition of graviton itself. For 
example, in addition to the two degrees of freedom associated to massless fields, there is also 
a scalar degree of freedom associated to the non-vanishing trace of the solution. 

On the other hand, when considering gravitational waves, since their sources are at enor¬ 
mous distance from Earth, the situation is somewhat different. To begin with, let us remark 
that, although plane waves have constant amplitude, the spherical symmetry of physical gravi¬ 
tational waves implies that their amplitudes fall off with the distance from the source. The am¬ 
plitude of the first-order solution, for example, scales with distance z according to A^^^u ~ I/ 2 ;. 
Considering that second-order effects are proportional to A^^^A^^^p, it is usually argued that 
second-order gravitational waves fall off with distance as 1 j ^, and are for this reason neglectful 
14l |. This is the case of the A-part of the second-order solution (I32p . whose amplitude falls off 
as 

However, due to an additional linear dependence on the distance, the amplitude of the H-part 
of that solution falls off as 

~ 1 / 2 :. (79) 

Contrary to the usual belief, therefore, second-order effects may be relevant even at large 
distances from the source. Neglecting the faster decaying A-part, the dominant solution at 
large distances is 

(80) 


hP 


"(2)1 


A2koXP 


From the second-order solution 
physical components are 


we see that this wave is transverse and traceless. Its two 


70 — 

"'( 2)0 ~ "'( 2)2 


and 


h'J — —h 
"( 2)2 “ "( 2)2 


(81) 


For astrophysical and cosmological applications, therefore, a kind of screening mechanism 
eliminates the traceful part of the wave. 

Finally, as a byproduct, we have studied the Kretschmann and Pontryaguin invariants. It 
is well-known that both invariants vanish identically for first-order gravitational waves. By 
using the higher-order solutions obtained in this paper, we have shown that these invariants 
vanish actually for gravitational waves of any order. 


11 




Acknowledgments 

The authors would like to thank FAPESP, CNPq and CAPES for financial support. 


References 

[1] See, for example, H. Bondi, Gravitational waves in general relativity, Current Contents 
30, 16 (1990), and references therein. Available at: 

http://www.garfield.library.upenn.edu/classicsl990/Al990DN22600001.pdf 

[2] M. Maggiore, Gravitational Waves (Oxford University Press, Oxford, 2008). 

[3] L. Blanchet, Gravitational Radiation from Post-Newtonian Sources and Inspiralling Com¬ 
pact Binaries, Liv. Rev. Rel. 5, 3 (2002). 

[4] D. Christodoulou and B. G. Schmidt, Convergent and asymptotic iteration methods in 
general-relativity, Commun. Math. Phys. 68, 275 (1979). 

[5] R. Aldrovandi, J. G. Pereira, R. da Rocha and K. H. Vu, Nonlinear gravitational waves: 
their form and effects, Int. J. Theor. Phys. 49, (2010); arXiv:0809.2911, 

[6] J. Ehlers and W. Kundt, in Gravitation: An Introduction to Current Research, ed. by E. 
Witten (Wiley, New York, 1962). 

[7] S. Hawking and G. Ellis, The Large Scale Structure of Space-Time (Cambridge University 
Press, Cambridge, 1973), page 260. 

[8] H. Stephani, General Relativity (Cambridge University Press, Cambridge, 1982), page 
164. 

[9] H.-J. Schmidt, Why do all the curvature invariants of a gravitational wave vanish? in New 
frontiers in gravitation, ed. by G. Sardanashvily (Hadronic Press, Palm Harbor, 1996); 
arXiv:gr-qc/9404037. 

[10] See, for example, R. K. Dodd, J. C. Eilbeck, J. D. Gibbon and H. C. Morris, Solitons and 
Nonlinear Wave Equations (Academic Press, London, 1982). 

[11] T. W. B. Kibble, The quantum theory of gravitation, in High-Energy Physics and Elemen¬ 
tary Particles (IAEA, Vienna, 1965), page 885. 

[12] The ALPHA Collaboration & A. E. Charman, Nature Communications 4, 1785 (2013). 

[13] R. Aldrovandi, J. G. Pereira and K. H. Vu, The nonlinear essence of gravitational waves. 
Found. Phys. 37, 1503 (2007); arXiv:0709.1603, 

[14] See, for example, B. Schutz, Gravity from the ground up (Cambridge University Press, 
Cambridge, 2003). 


12 



